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Abstract 

I discuss a recent application of homotopy perturbation method to a heat transfer 
problem. I show that the authors make infinity equal five and analyze the conse- 
quences of that magic. 

There has recently been great interest in the application of several approximate 
procedures, like the homotopy perturbation method (HPM), the Adomian de- 
composition method (ADM), and the variation iteration method (VIM), to 
a variety of linear and nonlinear problems of interest in theoretical physics 
[1-15]. From now on I will refer to those variation and perturbation approaches 
as VAPA. In a series of papers I have shown that most of the VAPA results 
are useless, nonsensical, and worthless [16-20]. In many of those papers the 
authors try to solve nonlinear problems by means of elaborated VAPA im- 
plementations that merely produce the Taylor expansion of the solutions. Of 
course, such approximate expressions do not give the overall picture of the 
dynamics, and the authors are merely content with a description of the ini- 
tial stages of the evolution which do not tell us anything relevant about the 
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process. Other authors solve the Schrodinger equation and obtain trivial un- 
physical solutions that are not square integrable. As an example I mention two 
great feats of the VAPA users: the expansion of exponential functions of the 
form e*"* [15] and a prey-predator model that predicts a negative population 
of rabbits [14] (see also my comments [18,19]). 

However, my criticisms have not been welcome because they lack "the qualities 
of significant timeliness and novelty that we are seeking in this journal" and 
for that reason they remain unpublished outside arXiv. 

The purpose of this article is the analysis of a recently published paper that 
certainly meets the criterion of timeliness and novelty sought in that journal. 
Esmaeilpour and Ganji [5] applied homotopy perturbation method (HPM) to 
the solution of the problem of forced heat convection over an horizontal flat 
plate. After some algebraic manipulation of the Navier-Stokes equations they 
obtained two coupled nonlinear differential equations: [5] 

r{v) + lmf"{v)-o 

with the boundary conditions 

/(0) = /'(0)=0, /'(oo) = l 

^(0) = 1, ^(oo) = (2) 
The HPM yields series of the form 

oo oo 

/-E/.p''^ = E^y (3) 

j=0 j=0 

where the dummy perturbation parameter p is set equal to unity at the end 
of the calculation. Esmaeilpour and Ganji [5] choose the boundary conditions 
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/,(0) = 0, /;(0) = 0, /j(oo)=(5,o 

%(0) = 5,o, ^,(oo)=0 (4) 



Surprisingly, the perturbation corrections fj{ri) and Oj{r]) are polynomial func- 
tions of 1] [5] which cannot satisfy the boundary conditions at infinity (jl]) al- 
though the authors appear to state otherwise [5]. In fact, the approximate 
function 

, , 1348969 2 4867 . 451 . 1 „ 

fupMiri) = r]^ ri^ H ri^ r?" (5) 

■> ^ 7741440 ' 10752000 ' 322560000 ' 532224000 ' ^ ^ 

corrected to third order (j < 3) does not satisfy the boundary conditions ([2]). 
However, the figures shown by Esmaeilpour and Ganji [5] exhibit a reasonable 
agreement between the exact and approximate solutions for < ?7 < 5. 

When VAPA does not fit the problem the users make the problem fit VAPA. 
In this case Esmaeilpour and Ganji [5] do some kind of magic and make infin- 
ity equal five. Consequently, their approximate solutions satisfy the following 
boundary conditions 

/,(0) = 0, /;(0) = 0,/j(5) = 5,o 

9^(0) = 6,0, 0^(5) = (6) 

Unfortunately, the authors forgot to say how they did this miracle. Since I am 
not that smart and still think that there is something else beyond that shrunk 
infinity I produced Fig. [1] that shows the actual behaviour of f'HPMiv) ^ 
wider interval. 

When solving the differential equation for / one has to determine the value of 
/"(O) that is consistent with the boundary condition at infinity. Esmaeilpour 
and Ganji [5] do not discuss the calculation of this unknown parameter al- 
though they obtained the numerical solution by a standard software. Our 
straightforward approximate calculation based on trial and error suggests that 
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/"(O) ^ 0.3320574 and the HPM function ([5]) yields /^Piv/(0) = 0.349. The dis- 
crepancy is probably due to the fact that I have not been able to enter the 
shrunk infinity discovered by the authors. I suppose that for this very rea- 
son my contribution cannot be considered to carry the qualities of significant 
timeliness and novelty. 

It is my opinion that VAPA have produced one of the greatest concentrations 
of bad papers I have ever seen. If the reader proves me wrong I will certainly 
apologize. 
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Fig. 1. Numerical (dashed line) and HPM (solid line) values of /'(r/) 
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